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Abstract 

We address the question of how hght can be the hghtest super- 
symmetric particle neutrahno to be a rehable cold dark matter (CDM) 
particle candidate. To this end we have performed a combined analysis 
of the parameter space of the Minimal Supersymmetric Standard Model 
(MSSM) taking into account cosmological and accelerator constraints 
including those from the radiative 6 — 57 decay. Appropriate grand 
unification (GUT) scenarios were considered. 

We have found that the relaxation of gaugino mass unification is 
sufficient to obtain a phenomenologically and cosmologically viable so- 
lution of the MSSM with the neutralino as light as 3 GeV. 

We have found good prospects for direct detection of these su- 
perlight CDM neutralinos via elastic scattering off various nuclei in 
the forthcoming experiments with low-threshold DM detectors. 

In a certain sense, these experiments can probe the gaugino mass 
unification giving constraints on the possible GUT scenarios within the 
MSSM. 



1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) is a leading candidate for a 
low energy theory consistent with the grand unification (GUT) idea. The gauge coupling 
constants precisely measured at LEP make unification in the Standard Model (SM) rather 
problematic while in the MSSM it occurs naturally with an excellent precision. 

The MSSM, supplied a priori with a complete set of the grand unification conditions, pos- 
sesses a remarkable predictive power. The complete set of the GUT conditions includes gauge 
coupling constants unification as well as unification of the "soft" supersymmetry (SUSY) 
breaking parameters at the same GUT scale Mx ~ 10^^ GeV. Instead of this uhimate 
GUT scenario one can consider less restrictive particular GUT scenarios relaxing some of 
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the GUT conditions. We have no yet firm theoretical arguments in favor of one of these 
scenario. Analyzing the prospects for discovering SUSY in various experiments it is more 
attractive to adopt a phenomenological low-energy approach and disregard certain GUT 
conditions. Following these arguments we will discuss several GUT scenarios in the present 
paper. 

Another advantage of the MSSM is the prediction of a stable lightest supersymmetric 
particle (LSP) - neutralino (x)- Now the neutralino is the best known cold dark matter 
(CDM) particle candidate 0. 

There is a well-known lower limit on the neutralino mass > 18.4 GeV p, This 
result is strongly connected to the unification scenario with the universal gaugino masses 
Mj(Mx) = Tni/2 at the GUT scale Mx- The renormalization group evolution of Mj, starting 
with the same mi/2 value, leads to tight correlation between the neutralino m^. (with i = 
1 — 4), chargino m^^ ^ gluino nig masses (see section 2). As a result, direct and indirect 
SUSY searches at Fermilab and at LEP Q strongly disfavor the neutralinos lighter than 
18.4 GeV within the universal gaugino mass scenario 0]. As discussed in the following, 
the non-universal gaugino mass scenario with non-equal gaugino masses at the GUT scale 
allows essentially lighter neutralinos. In a certain sense, direct searches for the "superlight" 
neutralinos in the mass range < 18.4 GeV could be a test of the gaugino mass unification. 

In the present paper we will consider the discovering potential of DM experiments search- 
ing for the superlight DM neutralinos via neutralino-nucleus elastic scattering. Since the 
nuclear recoil energy in this collision is E^. ~ lO^^m^, a DM detector should have the low 
threshold ~ few KeV to be sensitive to the DM neutralinos as light as ~ few GeV. 

There are several projects of DM experiments with low-threshold detectors which are 
able to probe this mass region ||^. Some of them are expected to run in the near future. 
These experiments will use either a new generation of cryogenic calorimeters [j^, |^, ^ or 
Germanium detectors of special configuration [1^. 



We will argue that these experiments have good prospects for verification of the gaugino 
mass unification at the GUT scale. The basic reason follows from the results of our analysis 
of the MSSM parameter space within the non-universal gaugino mass unification scenarios. 
In these scenarios we have found superlight neutralinos with masses as small as 3 GeV which 
produce cosmologically viable relic density 0.025 < ^^^^ < 1 and a substantial total event 
rate i? ~ 1 event/kg/day of elastic scattering from nuclei in a DM detector. These values of 
R are within the expected sensitivity of the above-mentioned low-threshold DM detectors. 
Therefore, the superlight DM neutralinos appearing in the non-universal gaugino mass GUT 
scenarios can be observed with these set-ups. Negative results of DM neutralino searches in 
the mass region 3 GeV< < 18.4 GeV would discriminate these scenarios. 

The paper is organized as follows. In section 2 we specify the MSSM and give formulas 
used in the subsequent sections. In section 3 we summarize the experimental and cosmologi- 
cal inputs for our analysis. Section 4 is devoted to calculation of the event rate of the elastic 
neutralino-nucleus scattering. In section 5 we discuss the results of our numerical analysis 
and section 6 contains the conclusion. 
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2 The Minimal Supersymmetric Standard Model 



The MSSM is completely specified by the standard SU(3)xSU(2)xU(l) gauge couplings as 
well as by the low-energy superpotential and "soft" SUSY breaking terms [|I|. The most 
general gauge-invariant form of the R-parity conserving superpotential is 

W = hEL'E'HUij + hoQWHieij + huQ'U'Hieij + fxHlHieij (1) 

(ei2 = +!)• The following notations are used for the quark (5(3,2,1/6), D'^{3, 1,1/3), 
[/"(3, 1, -2/3), lepton L{1, 2, -1/2), E%1, 1, 1) and Higgs Hi{l, 2, -1/2), 1/2(1, 2, 1/2) chiral 
superfields with the SU(3)cXSU(2)iXU(l)y assignment given in brackets. Yukawa coupling 
constants hE,D,u are matrices in the generation space, non-diagonal in the general case. For 
simplicity we suppressed generation indices. 

In general, the "soft" SUSY breaking terms are given by ||TT| : 



^SB = -lY.MAXAXA-mj,^\H,\^-mj,jH2f-ml\Q\^-ml\Df-ml\Uf 



^'a 

m||L|2 - m||E^|2 - {hEAEL^E'Hleij + hoADQWHle 



+ huAuQ'il^Hiei, + h.c) - {B^iH{Hiei, + h.c) (2) 

As usual, M^^2,i are the masses of the SU(3) xU(2) xU(l) gauginos g, W, B and rrii are the 
masses of scalar fields. Al, Ad, Ajj and B are trilinear and bilinear couplings. 

Observable quantities can be calculated in terms of the gauge and the Yukawa coupling 
constants as well as the soft SUSY breaking parameters and the Higgs mass parameter 
introduced in Eqs. (|l])-(||). 

Under the renormalization they depend on the energy scale Q according to the renormal- 
ization group equations (RGE). 

It is a common practice to implement the GUT conditions at the GUT scale Mx- It 
allows one to reduce the number of free parameters of the MSSM. A complete set of GUT 
conditions is: 

mi{Mx) = = mH,,,{Mx) = m^{Mx) = niQiMx) = m^{Mx) = mf,{Mx) = mo, (3) 
AuiMx) = AoiMx) = ALiMx) = Ao, 

Mi{Mx) = mi/2, (4) 

, , ^ , , 5g''^ gl 

ai{Mx) = acuT, where ai = - — , 02 = = — , (5) 

3 Att Air Air 

g' ,g and gg are the U(l), SU(2) and SU(3) gauge coupling constants. 

The above GUT conditions put very stringent constraints on the weak scale particle spec- 
trum and couplings. In this scenario the neutralino is heavier than 18.4 GeV, as mentioned 
in the introduction. Due to specific correlations between the MSSM parameters at the weak 
scale, stemming from Eqs. (^-(1), the detection rate of the DM neutralinos is very small, 
being typically beyond the realistic abilities of the present and future DM detectors. 
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However, at present we have no strong motivation to impose a complete set of the GUT 
conditions in Eqs. d)-® on the MSSM. 

Therefore, in our analysis we use relaxed versions of the above discussed GUT conditions 
and consider two scenarios with the non- universal gaugino masses M3 = M2 7^ Mi. 

(a) with the universal 

rriH^^iMx) = mj{Mx) = mo, and 

(b) non-universal 

rriH, (Mx) = niH, (Mx) ^ m^. (6) 

GUT scale scalar masses. Deviation from the GUT scale universality of the mass parameters 
in the scalar sector was previously considered in Refs. 

Accepting the GUT conditions, we end up with the following free MSSM parameters: 
the common gauge coupling acuT] the matrices of the Yukawa couplings hf' , where i = 
E,U,D; the soft supersymmetry breaking parameters mo, m-1/2, Aq, B, the Higgs field 
mixing parameter fi and an additional parameter of the Higgs sector rriA being the mass of 
the CP-odd neutral Higgs boson. Since the masses of the third generation are much larger 
than the masses of the first two ones, we consider only the Yukawa coupling of the third 
generation and drop the indices a, b. 

Additional constraints follow from the minimization conditions of the scalar Higgs poten- 
tial. Under these conditions the bilinear coupling B can be replaced in the given list of free 
parameters by the ratio tan/3 = V2/V1 of the vacuum expectation values of the two Higgs 
doublets. 

We calculate the weak-scale parameters in Eqs.(|T])-(0) in terms of the above-listed free 



parameters on the basis of 2-loop RGEs following the iteration algorithm developed in ||14 
The Higgs potential V including the one-loop corrections AV can be written as: 



2 '2 



with AV = ^5:(-l)2^>(2J, + l)amf 
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where the sum is taken over all possible particles with the spin Jj and with the color degrees 
of freedom Ci. The mass parameters of the potential are introduced in the usual way as 

<2 = <,,+/^^ ml = B^^, (8) 

They are running parameters mi{Q) with the Q-scale dependence determined by the RGE. 
The 1-loop potential (|^) itself is Q-independent up to, irrelevant for the symmetry break- 
ing, field-independent term depending on Q. At the minimum of this potential the neutral 
components of the Higgs field acquire non-zero vacuum expectation values {H^^} — "^1,2 
triggering the electroweak symmetry breaking with g'^{v1 + v^) = 2M^. 
The minimization conditions read 



2ml = 2m3tan/3-M|cos2/3-2Si, (9) 
2ml = 2m2cot/3 + M|cos2/5-2S2, (10) 
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where 



dAV 



, with iIji^2 = Re-ffi 2, are the one-loop corrections [|T^: 



327r' 



\2Ji 



(2.7, + 1) 



(11) 



As a remnant of two Higgs doublets Hi^2 after the electroweak symmetry breaking there 
occur five physical Higgs particles: a CP-odd neutral Higgs boson A, CP-even neutral 
Higgs bosons H, h and a pair of charged Higgses . Their masses , fj , mj^± can 
be calculated including all 1-loop corrections as second derivatives of the Higgs potential in 
Eq. d^) with respect to the corresponding fields evaluated at the minimum |17 . 



The neutralino mass matrix written in the basis {B, W^, H^, H2) has the form: 
/ 

= 



Ml —Mz cos P siuw Mz sin /? sinvi^ \ 

M2 Mz cos /? cosvy —Mz sin [3 cos^y 

Mz cos 13 sm\Y Mz cos (3 cosw — /x 

V Mz sin (3 sm.\Y —Mz sin {3 cosw — /i 



The universal gaugino mass unification scenarios with the GUT condition in Eq. 
the relation 

5 

Ml = - tan^ OwM2. 
3 



(12) 



imply 



(13) 



Diagonalizing the mass matrix by virtue of the orthogonal matrix M one can obtain four 
physical neutralinos Xi with the field content 



(14) 



and with masses m^. being eigenvalues of the mass matrix (|T2[). We denote the lightest 
neutralino Xi by X- our analysis x is the lightest SUSY particle (LSP). 
The chargino mass term is 



+ h.c. 



with the mass matrix 



M2 \/2Mw sin (3 

\/2Mw cos (3 fi 



which can be diagonalized by the transformation 

X- = UaW- + U,2H-, x^ = VaW+ + Vi2H' 
with U*J^^±V'^ = diag(M-±, M-±), where the chargino masses are 



(15) 



(16) 



(17) 



^1,2 



M^ + fi^ + 2M^T 



(M| - /i2)2 + 4M^ cos2 2(3 + AM^{Mi + /i^ + 2M2/i sin 2(3) 
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It is seen from Eqs. (|T^) and (0) that in the universal gaugino mass scenarios, leading to 
relation (0), the neutralino and chargino sectors of the MSSM are strongly correlated since 
they are described by the same set of three parameters M2,/i, tan/3. Relaxing condition 
(|13|) we get one independent parameter in the neutralino sector, Mi. In this case, taking 
place within the non-universal gaugino mass GUT scenarios, the neutralino mass can be 
appreciably smaller than 18.4 GeV. The latter is the lower neutralino mass bound following 
from the experimental data in the case of the universal gaugino mass. 

The mass matrices for the 3-generation sfermions t, b and f in the fi — fn basis are: 

m| + mf + \{^M^ - Ml) cos 2(3 mt(A - /i cot (3) \ 

mt{At- II cot 13) m| + m2-|(M^-M|)cos2/5 /' 
m| + - |(2M^ + M|) COS2/5 /^^(Ab - /itan/3) \ 

mb(v4fe - /itan/3) m|, + + |(M^ - M|) cos 2/3 /' 

m| + m2 - |(2M^-M|)cos2/3 m^(A^ - /itan/3) \ 

mr{Ar - ^itaIl|3) m| + + (M^ - M|) cos 2/5 )' 

For simplicity, in the sfermion mass matrices we ignored non-diagonality in the generation 
space, which is important only for the b —>■ s'j -decay. 



Ml = 
Ml = 
Ml = 



3 The Constrained MSSM parameter space 

In this section we summarize the theoretical and experimental constraints used in our anal- 
ysis. 

The solution of the gauge coupling constants unification (see Eq. (|^)) allows us to define 
the unification scale Mx- Our numerical procedure is based on the 2- loop RGBs. We use 
the world averaged values of the gauge couplings at the energy obtained from a fit to the 



LEP data [IS], Mw [H and rrit |19, E 



a~\Mz) 
sin^ Ojjs 

"3 



128.0 ±0.1, 
0.2319 ±0.0004, 
0.125 ±0.005. 



(18) 
(19) 
(20) 



The value of the fine structure constant a ^{Mz) was updated from ||2T| by using new data 
on the hadronic vacuum polarization |2^. The standard relations are implied: 



a 



a 



3 cos^ 6*^1/ 



012 



sin^^ 



w 



(21) 



SUSY particles have not been found so far and from the searches at LEP one knows that the 
lower limit on the charged sleptons is half the Z° mass (45 GeV) and the Higgs mass 
has to be above 60 GeV p3|, while the sneutrinos have to be above 41 GeV 0]. For 
the charginos the preliminary lower limit of 65 GeV was obtained from the LEP 140 GeV 
run P]. The above mass limits are incorporated in our analysis. 
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Radiative corrections trigger spontaneous symmetry breaking in the electroweak sector. 
In this case the Higgs potential has its minimum for the non-zero vacuum expectation values 
of the Higgs fields if°2- Solving Mz from Eqs. (p|)-(p^ yields: 

Mz^ _ m^ + Si-K + S2)tan^/? 
2 ~ tan2/3-l ' ^ ^ 

where Si and E2 are defined in Eq. ([111) . This is an important constraint which relates the 
true vacuum to the physical Z-boson mass Mz = 91.187 ± 0.007 GeV. 

Another stringent constraint is imposed by the branching ratio BR{b — 57) , measured 
by the CLEO collaboration ^ to be BR{b 57) = (2.32 ± 0.67) x 10"^ 

In the MSSM this flavor changing neutral current (FCNC) process receives contributions 
from H"^ —t,x^ — t and g — q loops in addition to the SM W — t loop. The x ~ t loops, 



which are expected to be much smaller, have been neglected The g — q loops are 



proportional to tan /3 . It was found that this contribution should be small, even in 
the case of large tan j3 and therefore it can be neglected. The chargino contribution, which 
becomes large for large tan (3 and small chargino masses, depends sensitively on the splitting 
of the two stop mass eigenstates ti 2- 



Within the MSSM the following ratio has been calculated |26|: 



BR{b .7) _ \V:M' ^QcnGa [v''^''A, + |(r^W23 _ ^16/23)^^ + c 



2 



BR{h^ceu) TT I{x,b)[l-{2/?,Ti)as{mb)f{x,b)] ' 

where 

0.175, /(xefe) = 0.4847, r/ = a,(MH/)/a,(mb), /(x^f.) = 2.41. 

Here f{xcb) represents corrections from leading order QCD to the known semileptonic b —>■ 
ceu decay rate; I{xcb) is a phase space factor; Xcb = rric/mf, = 0.316; K^^q describes the 
next-to- leading-order QCD corrections p8[. A^^g are the coefficients of the effective operators 
for 6s-7 and bs-g interactions respectively; C describes mixing with the four-quark operators. 



The ratio of CKM matrix elements '^'^^fJ" = 0.95 was taken from |29^ 



Assuming that the neutralinos form a dominant part of the DM in the universe one 
obtains a cosmological constraint on the neutralino relic density. 

The present lifetime of the universe is at least 10^° years, which implies an upper limit 
on the expansion rate and correspondingly on the total relic abundance. Assuming /ig > 0.4 
one finds that the contribution of each relic particle species x has to obey pO|: 



n^hl < 1, (23) 

where the relic density parameter = p^/ is the ratio of the relic neutralino mass density 
p^ to the critical one pc = 1.88 ■ lO^^^hgg-cm^'^. 

We calculate ^x^o following the standard procedure on the basis of the approximate 
formula 

,2 _ o 10 in-ii V ^7 ,,1/2 / GeV-2 




n^h' = 2.13xlO-^M^ xAT^^M • , . (24) 

° ' ^'^^T^ol F \axF + bxl/2J ^ ' 
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Here is the present day photon temperature, T^jT^ is the reheating factor, xp = Tpjra^ ^ 
1/20, Tp is the neutrahno freeze-out temperature, and Np is the total number of degrees of 
freedom at Tp- The coefficients a, b are determined from the non-relativistic expansion 



< CTannyhV >^ a + bx 



(25) 



of the thermally averaged cross section of neutralino annihilation. We adopt an approximate 
treatment ignoring complications, which appear when the expansion ( ^51) fails [0. In our 
analysis all possible channels of the x~X annihilation are taken into account. The complete 
list of the formulas, which we used, for the coefficients a, b and numerical values for the other 



parameters in Eqs. (^) and ( pSj ) can be found in ^ ^ 

Since the neutralinos are mixtures of gauginos and higgsinos, the annihilation can occur 
via s-channel exchange of the Z^, Higgs bosons and t-channel exchange of a scalar particle, 
like a selectron ||33|. Therefore, the cosmological constraint in Eq. (p3D substantially restricts 
the MSSM parameter space, as discussed by many groups p|, ^ ^ 

In the analysis we ignore possible rescaling of the local neutralino density p which may 
occur in the region of the MSSM parameter space where flyh'^ < 0.025 pBl p^. At lower relic 



densities DM neutralinos cannot saturate even galactic halos in the universe and the presence 
of additional DM components should be taken into account. One may assume that it can be 
done by virtue of the above-mentioned rescaling ansatz. Let us note that the halo density 
is a very uncertain quantity. Its actual value can be one order of magnitude smaller (or 
larger) than the quoted value 0.025 |38|. The SUSY solution of the DM problem with such 
low neutralino density becomes questionable. Therefore, we simply skip the corresponding 
domains of the MSSM parameter space as cosmologically uninteresting. 

Thus, we assume neutralinos to be a dominant component of the DM halo of our galaxy 
with a density = 0.3 GeV-cm~^ in the solar vicinity. 



4 Neutralino-Nucleus Elastic Scattering 

A dark matter event is elastic scattering of a DM neutralino from a target nucleus producing 
a nuclear recoil which can be detected by a detector. The corresponding event rate depends 
on the distribution of the DM neutralinos in the solar vicinity and the cross section a^iix^) 
of neutralino- nucleus elastic scattering. In order to calculate aei{x^) one should specify 
neutralino-quark interactions. The relevant low-energy effective Lagrangian can be written 
in a general form as 

Leff = T.(A- X7m75X ■ gT'^Tsg +^-C,-xX-qq) + O (^] , (26) 

where terms with vector and pseudoscalar quark currents are omitted being negligible in the 
case of non-relativistic DM neutralinos with typical velocities ~ 10~^c. 

In the Lagrangian ( PU] ) we also neglect terms which appear in supersymmetric models at 
the order of and higher, where rriq is the mass of the scalar superpartner q of the quark 
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q. These terms, as pointed out in Ref . , are potentially important in the spin-independent 
neutralino-nucleon scattering, especially in the domains of the MSSM parameter space where 
rriq is close to the neutralino mass m^. Below we adopt the approximate treatment of these 



terms proposed in which allows "effectively" absorbing them into the coefficients Cg in 
a wide region of the SUSY model parameter space. Our formulas for the coefficients Ag and 
Cq of the effective Lagrangian take into account squark mixing — qn. and the contribution 
of both CP-even Higgs bosons /i, H: 

^3 





1^ Id 




L 2 






m?i - 









cos^ 9g + sin^ 9q 




2vsin 9q 



+ COS^ Qq 



+ 



+ mq) 



Mw Pg sin 2^, T3(A/'i2 - taxiOwUu 



m'^qi - (m^ + mqY' 



[m^ + mq) 



(27) 



Here 



2 <? 2 <? 

4 ^mt m 



+ Pa 



ml 

( COS^ Ba 



H 



sin^ On 



COS^ da 



sin^ Oa 



')qL 



m 



+ sin 



X 

Fh 

hq 
Hg 

(t^qR 
P. 



<?1 



- (m^ + mqY 



1 



-P 



w 



m„ 



V (PqR) 



- i^x + ^qY - i^x + ^q) 



(7Vi2 — Ml tan 9w) (A^u cos an + Ms sin an), 
(A/i2 — A/ii tan 9w) (A/i4 sin an — A/13 cos an), 



, 1 , cos 
2 smp 



1 ^ sin an 

2 cosp 



1 sina/f 1 cosa// 

2 smp 2 cosp 
AAi2T3+Arn(g-T3)tan^^, 

tan6'vK Q A/11, 

(i + T3)^ + (i-T3)% 

2 smp 2 cosp 



(28) 
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The above formulas coincide with the relevant formulas in ||32| neglecting the terms ~ l/m| 
and higher. These terms are taken into account "effectively" by introducing an "effective" 
stop quark t propagator. 

A general representation of the differential cross section of neutralino-nucleus scattering 
can be given in terms of three spin-dependent J^ij{q^) and one spin-independent J-'s{(f') form 
factors [^1 

da , n 



8Gf 



dq^ 



[v,q 



(29) 



where f is a projectile neutralino velocity and q is the momentum transferred to the nucleus. 
The last term corresponding to the spin-independent scalar interaction gains coherent en- 
hancement {A is the atomic weight of a nucleus in the reaction). The coefficients ao,i, Cq 
do not depend on the nuclear structure and relate to the parameters Aq,Cq of the effective 
Lagrangian (p6|) as well as to the parameters Ag, / characterizing the nucleon structure. 
One has the relationships 



ao 
Co 



{Au + Ad){Au 
{Au-Ad){Au 



Ad) + 2AsAs 
Ad), 



f^^^^^ + fsC. + hl-fs-f)iC. 
TUu + rud 27 



Cb + Ct 



(30) 



Here Ag^^"-' are the fractions of the proton (neutron) spin carried by the quark q. The 
standard definition is 

< pin)\qY7,q\pin) >= 2S;^^^Aq^^-\ (31) 

where S^^^-^ = (0, Sp{n)) is the 4-spin of the nucleon. The parameters Ag^*^"^ can be extracted 
from the data on polarized nucleon structure functions ^ and the hyperon semileptonic 
decay data p2| . 

We use Ag values extracted both from the EMC |^ and SMC |^T| data. The other 



nuclear structure parameters fs and / in formula (^) are defined as follows: 

<p{n)\{mu + md){uu + dd)\p{n)> = 2/Mp(„)^^, 
< p{n)\msss\p{n) > = fsMp(^n)^^. 



(32) 



The values extracted from the data under certain theoretical assumptions are P3| : 

/ = 0.05 and /, = 0.14. (33) 
The strange quark contribution fs is known to be uncertain to about a factor of 2. Therefore 



we take its value within the interval 0.07 < fs < 0.3 ^ 



The nuclear structure comes into play via the form factors ^ij(g^), ^^(g^) in Eq. (pQ]). 
The spin-independent form factor J^si^f) can be represented as the normalized Fourier trans- 
form of a spherical nuclear ground state density distribution p(r). We use the standard 
Woods-Saxon inspired distribution ESl. It leads to the form factor 



dhp{r)e"''^ 



qRo 



(34) 
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where Rq = {B? — bs^Y^"^ and s ~ 1 fm are the radius and the thickness of a spherical nuclear 
surface, ji is the spherical Bessel function of index 1. 

Spin-dependent form factors J-'ij{q^) are much more nuclear-model-dependent quantities. 
In the last few years noticeable progress in detailed nuclear-model calculations of these 
form factors has been achieved. For many nuclei of interest in the DM search they have 



been calculated within various models of nuclear structure Unfortunately, these 

calculations do not cover all isotopes which we are going to consider in the present paper. 
Therefore, we use a simple parameterization of the dependence of J^ij{q^) in the form of a 
Gaussian with the rms spin radius of the nucleus calculated in the harmonic well potential 



For our purposes this semi-empirical scheme is sufficient. 
An experimentally observable quantity is the differential event rate per unit mass of the 
target material 



dR 





/ 




Jv 



dvf{v)v—{v,Er), (35) 



dcr 

where is the nuclear recoil energy. The function f{v) is the velocity distribution of 
neutralinos in the earth's frame. In the galactic frame it is usually assumed to have an 
approximate Maxwellian form. Vmax = Vesc ~ 600 km/s and = 0.3 GeV-cm~^ are the 
escape velocity and the mass density of the relic neutralinos in the solar vicinity; Vmin = 
{MAEr/2Mla)^^^ with Ma and Mred being the mass of nucleus A and the reduced mass of 
the neutralino-nucleus system, respectively. Note that = 2MAEr. 

The differential event rate is the most appropriate quantity for comparing with the ob- 
served recoil spectrum. It allows one to take properly into account the spectral characteristics 
of a specific DM detector and to discriminate a background. However, discussing general 
problems and prospects of DM detection it is enough to consider the total event rate R 
integrated over the whole kinematical domain of recoil energy. Notice that this quantity is 
less sensitive to details of the nuclear structure than the differential event rate in Eq. (pS]). 
The q^ shape of the form factors J^ij{q^), J^si.'f) in Eq- ( p9|) may substantially change from 
one nuclear model to another. Integration over g^, as in the case of the total event rate i?, 
reduces this model dependence. 

The present paper aims at the general investigation of detectability of the cosmologically 
viable "superlight" DM neutralinos independently of a specific DM detector. Therefore, we 
may use the total event rate i? as a characteristic of the DM signal. 



5 Numerical Analysis and Discussion 

In our numerical analysis we randomly scanned the MSSM parameter space including all 
experimental and cosmological constraints discussed in section 3. Two GUT scenarios with 
the non-universal gaugino mass M3 = M2 7^ Mi (see Eqs. (§])-(§) and the comments below 
them) have been considered: (a) with the universal mHi2{Mx) = ^fi^x) = n^-o and (b) 
non-universal rriH^^Mx) = mH2{Mx) 7^ tiiq GUT scale scalar masses. 

The free parameters of both scenarios and intervals of their variation in our analysis are 
given in table I. 
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scenario 


Parameter 


Lower and Upper Bounds 


a,b 


tan/? 


1 


50 


a,b 


mo 


TeV 


1 TeV 


a,b 




-1 TeV 


1 TeV 


a,b 




-1 TeV 


1 TeV 


a,b 


M2 


TeV 


1 TeV 


a,b 


Ml 


-1 TeV 


1 TeV 


b 


rriA 


50 GeV 


1 TeV 



TABLE L The MSSM parameters within the GUT scenarios (a), (b) 
and intervals of their variations in the present numerical analysis. 

After the scan procedure had been finished we found the superlight neutralinos in the 
mass range 3 GeV < < 18.4 GeV in both GUT scenarios (a) and (b). The criterion we 
used to stop the running numerical code was stabilization of the lower bound. 

Recall that 18.4 GeV is the widely cited lower neutralino mass bound valid within the 
complete GUT scenario with relations d^)-® at the GUT scale. Scenario (a) implies a min- 
imal relaxation of the complete GUT scenario necessary to obtain the superlight neutralinos. 
Scenario (b) does not change the lower bound 3 GeV, but allows a much larger DM event 
rate R than scenario (a). We have calculated the total event rate for germanium (^^Ge), 
sapphire (AI2O3), fiuorine (^^F) and sodium iodide (Nal). 

The main results of our analysis are presented in Figs. 1-5 in the form of scatter plots. 
In Figs. 2-4 we cut off all points with R < 0.01 events/kg/day since they lie beyond the 
sensitivity of the present and the near-future DM detectors. 

Figure 1 shows the relic density, Q^h"^, produced by the superlight neutralinos and their 
gaugino fraction, Zg = A/'fi + J^u (see Eq. (p!4D). It is seen that the majority of the 
points are concentrated in the region of a large relic density, especially within the interval 
8 GeV < < 16 GeV. However, the neutralinos with a marginal mass of 3 GeV can 
also produce cosmologically interesting values of Q^h"^. Therefore we conclude that the 
superlight neutralinos can comprise a dominant part of the CDM. The field composition of 
these neutralinos characterized by the gaugino fraction Zg shows specific tendencies. A large 
domain of the MSSM parameter space (the volume is proportional to the number of points 
in the scatter plots) corresponds to the mixed Zg ~ 0.7 gaugino-higgsino states around the 
mass value 12 GeV. For larger neutralinos lose their higgsino component and become 
mostly gauginos. 

In Figures 2-3 we present the calculated event rate, R, for ^'^Ge, AI2O3, ^^F and Nal versus 
the neutralino mass in two GUT scenarios (a) and (b). These figures demonstrate that 
in both scenarios there are a lot of points within the reach of the near-future low-threshold 
DM detectors mentioned in the introduction. Their sensitivities are expected to be at the 
level of i? > 0.1 events/kg/day. In the scenario (b) points extend up to i? ~ 300 -7- 500 
events/kg/day for Ge, Nal and i? ~ 20 ^ 50 events/kg/day for ^^F, AI2O3. In scenario (a) 
these values are reduced by approximately a factor of 10. 

In unification scenario (b) with non-universal Higgs mass parameters (see Eq. (^) the 
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FIG. 1. The neutralino relic density Qh"^ and the gaugino fraction Zg = N"li + ■N'12 versus 



the neutrahno mass in GUT scenario fa). 



Higgs and sfermion masses are not strongly correlated. As discussed in |13|, this relaxation 



of the complete unification in the scalar sector makes it possible to avoid one of the most 
stringent theoretical limitations on the allowed values of the DM neutralino event rate. In 
this case the masses of the CP-even and charged Higgses mH,h, are calculated in terms 
of the CP-odd Higgs mass, uia, as an extra free parameter. 

An important question touches upon the fraction of the spin-dependent part Rsd of the 
total event rate R = Rsd + Rsi- From Eqs. (|27| ) -(|30| ) it is seen that measurement of Rsd 
and Rsi would give us complementary information about the MSSM parameters. In the 
previously investigated mass region > 18.4 GeV most of the experimentally interesting 
isotopes are more sensitive to the spin-independent part Rsi due to the coherent enhancement 
effect (see and references therein). The spin-dependent part Rsd originates from the 
neutralino-nucleus interaction via spin-spin coupling. It is not a coherent interaction with a 
whole nucleus since only a few nucleons contribute to the nuclear spin. 

In Figures 4 and 5 we give the scatter plots in i? - Rsd/ Rsi and Rsd/ Rsi - f^x planes for 
large event rate scenario (b). The following conclusions are valid in scenario (a) as well. One 
can see that in the mass region of the superlight neutralinos 3 GeV < < 20 GeV the 
coherent Rsi component of the event rate dominates for all isotopes in question, except the 
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FIG. 2. The expected DM neutralino event rate R versus the neutrahno mass for various 
isotopes of experimental interest. GUT scenario (a). All points with R < 0.01 events/kg/day 
are cut off. 



region of very small total event rates R. As seen from Fig. 5, the spin-dependent component 
Rsd dominates for the majority of the points in the mass range 8 GeV< < 14 GeV. It 
reflects the fact illustrated in Fig. 1 that in this mass region the neutralino can acquire the 
largest higgsino admixture 1 — f» 0.35. As a result, the Z-boson contribution to Rgd (via 
the coefficient Aq in Eq. (^71) ) is enhanced. 

We conclude this discussion with the following remark. Nuclear spin does not play an 
important role in the direct searches for the DM neutralinos in the mass region 3 GeV< < 
18.4 GeV. Previously the similar conclusion was obtained concerning the mass region 18.4 
GeV< |13[. Nevertheless, only a DM detector with a spin-non-zero target nuclei can 
provide us with information about Rgd and the corresponding MSSM parameters in Ag (see 
Eq. (g3)). 

6 Conclusion 

We have analyzed the MSSM parameter space taking into account cosmological and accel- 
erator constraints including those from the radiative b —>■ s'j decay. 

It is well known that the MSSM with the universal gaugino mass at the GUT scale 
disfavors neutralinos lighter than 18.4 GeV P, |^ if the known experimental data are taken 
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FIG. 3. The same as in Fig. 2 but for large event rate GUT scenario (b). 



into account. 

A central result of the present paper is the conclusion about the existence of a substantial 
domain of the MSSM parameter space corresponding to the superlight neutralinos in the mass 
range 3 GeV< < 18.4 GeV within the GUT scenarios with the non-universal gaugino 
mass. In this domain neutralinos are cosmologically viable and produce an event rate which 
is detectable in the near-future experiments with the low-threshold DM detectors. 



When our paper was being prepared, we found Ref. where prospects for the su- 
perlight DM neutralino detection are also considered but within a more phenomenological 
approach ignoring all GUT conditions. Our analysis shows that we need not disregard all 
GUT conditions in Eqs. to get a window for the superlight neutralinos. To this 

end it is necessary and sufficient to relax only gaugino mass unification condition. Since we 
know a generic root of the superlight neutralino we can conclude that in a certain sense the 
experiments searching for the DM neutralinos in the mass range > 18.4 GeV probe the 
gaugino mass unification at the GUT scale. 
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FIG. 4. Spin sensitivities of various isotopes of experimental interest. Given is the total event 
rate R in GUT scenario (b) versus the ratio Rsd/Rsi- Rsd and Rsi are the spin-dependent and 
spin-independent components oi R{R — Rsd + Rsi)- All points with R < 0.01 events/kg/day 
are cut off. 
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